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Abstract 

We study the dispersion relations of mesons and baryons built from Brillouin quarks on 
one Nj = 2 gauge ensemble provided by QCDSF. For quark masses up to the physical 
strange quark mass, there is hardly any improvement over the Wilson discretization, if 
either action is link-smeared and tree-level clover improved. For quark masses in the range 
of the physical charm quark mass, the Brillouin action still shows a perfect relativistic 
behavior, while the Wilson action induces severe cut-off effects. As an application we 
determine the masses of the 0°, f2+ and baryons on that ensemble. 



1 Introduction 

Wilson fermions offer an effective way of regulating the quark sector of QCD. Their conceptual 
simplicity entails a one-to-one correspondence between lattice and continuum flavor, a property 
which is particularly convenient for studying flavor physics in the Standard model or in one of 
its extensions. The main disadvantage is that they induce cut-off effects in physical observables 
which are both parametrically and numerically large, i.e. oc a, where a is the lattice spacing, 
and the prefactor is sizable (for a discussion see e.g. the recent review [I]). 

Two technical remedies which have proven useful in mitigating the discretization effects 
are clover improvement and link smearing. The first one changes the parametric behavior to 
anything between a n a and a 2 (a is the strong coupling constant) [2], depending on whether the 
coefficient csw in @ is adjusted in perturbation theory or non-perturbatively [3J. The second 
one concerns the links that enter the Dirac operator, in the covariant derivative and/or the 
clover term, and reduces the coefficient that multiplies the cut-off terms [4|j8]. It turns out that 
either idea greatly enhances the others effectiveness [9j; by combining a generic overall link 
smearing with an un-sophisticated improvement strategy (e.g. the tree-level choice csw — 1) 
the amount of chiral symmetry breaking becomes as small as am res 



In recent years such tree-level improved fat-link Wilson fermions have proved extremely 
successful, in particular in enabling simulations of QCD with Nf — 2 + 1 dynamical fermions 
directly at the physical mass point [IT] (for an overview of the physics results obtained with 
such studies see [l]). One can ask whether further action improvements would warrant the 



potential increase in CPU time needed to solve the Dirac equation Dx = b for a given right- 
hand side b. The goal of this paper is to investigate this question for the case of the "Brillouin 
fermion" proposed in [12]. From the definition of the Wilson and Brillouin actions 

D™(x,y) = E7,Vf(x,y)-^/A std (x,y) + ^-^E^A,, (1) 

D B «(x,y) = ETmVJT^I/) -hA bA (x,y) + y x>y -^^2a^F,J x>y (2) 

it is clear that the new fermion discretization is structurally close to the Wilson variety, except 
that the covariant derivative and the Laplacian are chosen to minimize the amount of rotational 
symmetry breaking (which led to an operator which is quite close to the one of Ref. 13 , in 
spite of the construction being based on rather different principles). 



The remainder of this article is organized as follows. In Sec. [2] we specify the set of Nf = 2 
gauge field configurations that we use to carry out our investigation, and we give further 
details of the link smearing and clover improvement (which we use both in the Wilson and in 
the Brillouin action, to compare like with like). Next, Sec. [3] contains the precise form of the 
Wuppertal smearing that we apply on both the source and the sink side of our propagators, and 
describes the procedure by which we tune the mass parameters k in ([II 12]) to the correct value for 
the light, strange and charm quark mass. Sec. [4] contains the central piece of our investigation, a 
comparison of the dispersion relation -E'(p) 2 as a function of the spatial momentum p 2 for mesons 
and baryons built from Wilson and Brillouin fermions. As a phenomenological application, we 
compare in Sec. [5] the mass of the Vt° c baryon that we find to experiment, and we give the masses 
of the hitherto unobserved states f2+ and on the ensemble considered. We summarize our 
findings in Sec.lHl and arrange details of the baryon interpolating fields in an appendix. 



2 Ensemble and valence action details 

The goal of our investigation is to compare the Wilson Q and Brillouin (|2]) fermion actions 
in the valence sector, with special emphasis on the dispersion relation E 2 = E 2 (p 2 ) for mesons 
and baryons composed of such quarks. We shall use a freely available set of dynamical gauge 
field configurations, i.e. with the effect of light sea quark loops included. We select an ensemble 
out of the Nf = 2 collection by QCDSF (Ml, namely the one with 



= 5.29, L/a = 40, T/a = 64, <f = 0.13632 (3) 



and aM* a = 0.1034(8) [14]. Thus M* ea L = 4.136, which bears the promise of small finite-size 
effects. QCDSF determines the scale of their j3 = 5.29 ensembles to be a -1 = 2.71(2)(7) |15|, 
tantamount to a = 0.0728(5)(19) fm, which implies M" ea ~ 280 MeV and L ~ 2.9 fm. 

Either action involves 3-fold APE smeared gauge links, where one smearing step in d — 4 
space-time dimensions is given by 

V^x) = Psu(3){(l-a)U„(x) + £ U v {x)U,{x + D)Ui{x+ft} (4) 



12 



and we use ct^ = 0.72. Here Psu(3) denotes the back-projection to SU(3) as described in 
and U is the unsmeared (original) gauge field. Note that these four-dimensionally smeared links 
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Figure 1: Left: 1 — (ReTr([/ spa t)/3) versus the number of 3D APE smearings for various values 
of a 3 D,' the perturbative bound is « 3 d < 2/3 [9]. Right: contour plot of (r 2 ) 1 / 2 of a Wuppertal 
smeared quark source; the asterisks refer to Fig.[2| Either plot is based on 5 configurations. 



enter both the (relevant) covariant derivative and the (irrelevant) clover term of the Wilson 
and Brillouin operators (for a summary of the options see [9], but this choice is the simplest 
and, as far as we can see, the most effective one). In addition, the clover coefficient is set to 



its tree- level value (csw — 1) for either action. In short, this is the same action as used in 12 
except that we now use three steps of APE smearing rather than one. 

3 Wuppertal smearing and quark mass tuning 

Wuppertal smearing amounts to the spreading of a vector g(x, t) - with non-trivial support on 
the time-slice t - within that time-slice by means of iV\y operations of the form [16] 

— > 1 , U — ? y.i (5) 

with spreading parameter <5w- The index i runs over the three spatial directions, and an implicit 
summation over y takes place. Here the spatial links Vi(x, t) are to be generated with the 
three-dimensional version of the APE smearing Q applied to the time-slice t. After the Dirac 
equation Dp = q has been solved with the broadened source q, an identical spreading is applied 
to the solution p (for each time-slice separately), i.e. we use smeared-smeared propagators with 
the same smearing on the source and on the sink side. 

In principle Wuppertal smearing has four parameters to adjust, namely (Nsr>, q;3d) of the 
three-dimensional APE smearing that the spatial links Vi(x, t) have undergone, and (Nw,Sw) 
in the recipe (g. For convenience we try to optimize the two pairs separately. 

We first study the behavior of the spatial plaquette of the lattices ^ under repeated 
applications of the three-dimensional version of the APE recipe Q. Our results are displayed 
in the left panel of Fig.[TJ In [9] there is the perturbative stability bound (d—l)/d in 

d space-time dimensions. Our results suggest that any a^" E < a^g 1 induces, asymptotically, 
a power-law fall-off of 1 — P3D, and thus allows to drive the 3D-plaquette arbitrarily small. 
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Figure 2: Effective mass of the pseudoscalar cc meson made from Brillouin fermions with 
k = 0.112429 for three different widths (cf. Fig.^ of the Wuppertal smeared sources and sinks. 



In addition, cn^o seems not too far from the perturbative prediction of 2/3 191. Hence, our 
recommendation is to apply a large number of 3D APE smearings on the spatial links that 
enter the Wuppertal spreading (|5]), e.g. (A^ 3 d,Q!3d) = (300,0.6) or (1000,0.6). 

To decide on the second pair (jVw,^w) m © we first consider the root-mean-square (rms) 
radius of the smeared source as a function of these parameters, the former being defined through 

2 = S r rV(r)g(r) 
E r ? t (r)g(r) 

Results from 5 configurations are shown in the right panel of Fig.[T} As expected, the contours 
look like hyperbolas, that is the rms radius of q is in the first place a function of the product 
iVw^w- Next, we consider the effective mass of a meson made from k = 0.112429 Brillouin 
fermions (which is in the vicinity of K^ n , see below) for (N w = 10, 30, 90, 5<w = 3), with results 
presented in Fig.[2j In this range, a higher iteration count appears not to decrease the coupling 
to excited states, but rather seems to induce more noise in the correlators. As a consequence 
we decide to stay with the conservative parameter set (iVw = 10, #w = 0.5). 

To compare like with like, we wish to compare the dispersion relation for mesons and baryons 
put together from Wilson and Brillouin fermions at a fixed value of the light, strange or charm 
quark mass. This is most conveniently done by first tuning the two ^-values to get common 
values of aM vl aM sS and aM c5 , respectively, for either discretization. For the mesons we use 
the PP correlators, and given their cosh-form it is advantageous to define the effective mass as 



aM eS (t) = - log ^ V ( ? ) 

2 vc(m) + Jc(t+iy -c(t/2) 2 < 



because this modification remedies the fall-off that the effective mass would show near the 
center of the box if the generic definition 

aM eS {t) = -log 1 



2 & VC(t + l) 
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Figure 3: Left: Tuning of K\ to achieve (aMp) 2 = 0.01069 with either action. Right: Similar 
tuning of k s to achieve (aM P ) 2 = 0.06404. 



would be used. For baryons with (exact or approximate) projection to a definite parity (see 
Sec. [4] below) we use the latter form, since these correlators do not show the cosh form. 

For the light quark we demand the pion mass to be the same (in lattice units) as in the 
sea, that is (aMp) 2 = 0.01069, given the information provided beneath ([3]). We solve the Dirac 
equation for a few K-values in the vicinity of the suspected target value, and interpolate them 
linearly to obtain the desired Ki, as shown in the left panel of Fig. [3} The results read 

Kf ri = 0.125249 , K^ il = 0.126146 . (9) 

For the strange quark we use the scale provided below ^ and the value M sS = 0.6858(7) GeV 
which follows via w2M|- — M 2 from the isospin-averaged and electromagnetically corrected 



masses M n = 134.8(3) MeV and M K = 494.2(5) MeV 17 . This gives the target value (aM P ) 2 = 
0.06404, and a similar interpolation procedure, shown in the right panel of Fig.|3j yields 

n® vi = 0.1251902 , n = 0.1243560 . (10) 

For the charm quark we proceed analogously to the strange case, except that we now use the 
value M cg = 2.9810(11) GeV from PDG This yields the target value (aM P ) 2 = 1.210, and 
with essentially the same kind of procedure we find 

nf l = 0.112336 , = 0.112513 . (11) 

A careful look at the left panel of Fig. [3] reveals that one of our trial Kr values for the 
Brillouin action happens to be rather light; for this point we find aM n = 0.067(5), tantamount 



to M w ~ 180 MeV. Following 11 19 we monitor the inverse iteration count of the solver 
(which is a proxy for the smallest eigenvalue of D^D) to make sure that we do not run into 
an "exceptional configuration" problem. The results of this monitoring are displayed in Fig.|4j 
Even for the lightest quark mass the distribution is roughly Gaussian, and the origin is 3 to 4 
standard deviations away from the median. This shows that even in this strongly non-unitary 
regime (where M™ 1 is about 100 MeV lighter than M^ ea ) the Brillouin operator can be safely 
inverted. We consider this an encouraging sign of the great stability of our action against 
fluctuations of the small eigenmodes, and think that this stability has the potential to render 
the Brillouin operator a cheap alternative to overlap or domain-wall fermions. 
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l/ n CG l/ n BiCGStal> 



Figure 4: Iteration count of the inversions with the trial K\ of Fig.[3] for the Brillouin action, 
using the CG (left) and the BiCGstab (right) algorithm for the lightest and all but the lightest 
masses, respectively. In each case 0(25) configurations are used. 



4 Meson and baryon dispersion relations 



With the tuned Ki, k s , k c of (|9j [lOj in hand, we are now in a position to study the dispersion 
relation E 2 = E 2 (p 2 ) for mesons and baryons composed of either Brillouin or Wilson fermions. 
To this end we consider two-point correlators of the form 

C M (t,p) = ^(J M (x,t)J A/ (O,0))e ipx (12) 

X 

C±(t,p)= TVi(l± T4 ) £<J B (x,t) J B (O,0))e ipx (13) 

X 

where M and B define the quantum numbers of the meson or baryon, respectively. For non-zero 



momentum the correct parity projector in (13) would read |(1 ± H74) 20,21 , but we stay 
with the simpler form, since we are always interested in the lower-mass parity partner (which 



spinor index, say u. and v. The projection to a definite spin can be done with 22 

1 



requires no projection at all). Finally, the sink and the source in (13) contain an uncontracted 
nd v. r . 



{pll 2 : 



3 r ■ .^(invPv-VvPivf) 
1 1 



1 



3 p 



3p< 

2 P»Pu + 



3p 



2 ijhlJPv ~ Pnlvi>) 



(14) 
(15) 



and similar expressions for P^l 21 22 as given in eqn. (9) of 

1 



22 



For p = they simplify to 



Si 



-Jilj 



(PI 



l/2x 
1 )ij 



,Hl3 



(16) 

(17) 



which are easy to implement and good enough for many purposes (we are always interested in 



the lowest-mass state). For details of (more advanced) spin projection see 23 . 
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Figure 5: Effective mass plots of the pseudoscalar (left) and vector (right) cc meson before 
(top) and after (bottom) averaging over the k- configurations that contribute to k 2 = 16. 



Having defined the correlation functions, we can now move on to the dispersion relations. 
We consider spatial momenta p = 2tt/L ■ k with < k 2 < 20. Thanks to the cubic symmetry 
among the spatial directions, in general several k-configurations contribute to a given k 2 . An 
effective mass plot before and after taking an average over the various contributions is shown in 
Fig.[5j for the pseudoscalar and vector meson, in the case of k 2 = 16. We see no big discrepancies 
before the average is taken, hence the averaging seems justified. 

Repeating this procedure for all k 2 yields the data presented in Fig. [6} We show the dis- 
persion relations for the pseudoscalar and the vector meson with quark content ss (top), sc 
(middle), cc (bottom). Evidently, for the heavier masses there is a significant difference be- 
tween the Wilson data (red circles) and those with the Brillouin discretization (blue squares). 
The full black line is not a fit, but the relativistic dispersion relation E 2 = p 2 + M 2 , starting 
from the first data-point (where the two actions were tuned to yield the same result for the ss 
and cc pseudoscalar states, but not for the remaining four states). This line shows that the 
discretization effects are induced by the Wilson action; within statistical errors the Brillouin 
data are free from such effects. 

Similarly, we can work out the dispersion relations for baryons. The data for the "decuplet- 
type" states fij ss , fi^J [which form a 20-plet under SU(4)], as well as for the "octet-type" 
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Figure 6: Dispersion relations for the pseudoscalar (left) and vector (right) meson with ss 
sc (middle), cc (bottom) quark content. The black line shows the relativistic E 2 = p 2 + 
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states fic SS , &>tcs [which form a 20'-plet under SU(4)] are shown in Fig.JTj To make it clear 
which states we consider, the interpolating fields of both the 20-plet and the 20'-plet are listed 
in the appendix. Again, the full black line is not a fit but the relativistic dispersion relation 
E 2 = p 2 + M 2 , and we stress that no further tuning of k- values was performed. Just like in 
the meson case, we see no significant difference in the regime of the physical strange quark 
mass (which explains why we refrain from looking at even lighter k values, such data would 
just be more noisy). However, with every strange quark that is replaced by a charm quark, the 
difference becomes more pronounced, up to the point where the dispersion relation of the 
is seriously distorted with Wilson fermions, but relativistically correct with Brillouin fermions. 
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Figure 8: Effective masses of the fl sss , fi ccc (decuplet-type states) and fl ssc , f2 scc (octet-type 
states) at p = for the Wilson (left) and Brillouin (right) action. 



combination 


Wilson 


Brillouin 






0.616(8) 


0.604(8) 


i^ccs 


f^css) /^sss 


0.617(6) 


0.612(6) 


i^ccc 


f^ccs) /^sss 


0.633(5) 


0.633(5) 


i^ccs 


^css) / ^css 
^ccs) /^ccs 


0.383(3) 
0.284(1) 


0.382(3) 
0.286(1) 



Table 1: Relative mass splittings between the various Q states, as determined with the Wilson 
(left) and the Brillouin (right) action. Here, the symbol fl xyz is meant as a shorthand for Mq x s . 



5 Masses of multiply charmed baryons 

As a byproduct of our investigation, and with the goal of spurring further improvements, we 
can quote the masses of the Q baryons considered in the previous section. Let us emphasize 
that this activity is based on the single ensemble ([3]). We will try to minimize the systematic 
effects on the numbers given below, but there is no way of reliably assessing the size of such 
systematic effects, given the data that we have. 

Let us begin by showing the effective mass plots of the four states fi~ ss , £^ s , ^tcsi ^tcc 
in Fig. [8] This is at p = 0, and we see no significant difference between the Wilson and the 
Brillouin action. Let us recall that the tuning was done in the meson sector, so this is a first 
non-trivial observation. To minimize the systematics it is usually a good idea to look at mass 
splittings, and to form dimensionless ratios. The former trick mitigates the effect of excited 
states contaminations and possible finite volume effects, the latter one reduces the sensitivity 
to the overall scale. The effective mass plots for all adjacent mass splittings, based on ratios of 
correlators like (C css )/ (C sss ), are shown in Fig.[9] Even at this level of zoom, the quality of the 
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Figure 9: Same as Fig.$ but for the splittings M Qca -M nssa and M Qcc -Ma caa and M Qcc -M Uct . 



data appears rather good, and we see only very mild differences between the Wilson and the 
Brillouin data. We list the relative mass splittings (normalized both with the lower one of the 
two states involved and with the Vl~ = Vl sss base state) in Tab. [I] 

Regarding phenomenological numbers, we should first mention that the experimental mass 
of the fi- (i.e. the J p = 3/2+, c = 0, s = 3 state in Tab.g is 1672.45(29) MeV, and the mass 
of the tt° c (i.e. the J p = 1/2+, c = 1, s = 2 state in Tab.g) is 2697.5(2.6) MeV [l8). Hence, with 
the ratios listed in Tab.[lj we can only compare the mass of the baryon to experiment, but 
not the masses of the f2+ (i.e. the J p = l/2 + , c = 2, s — 1 state in Tab.[2| and of the (i.e. 
the J p = 3/2 + , c = 3, s = state in Tab.[3|. We shall use the first three lines in Tab.[lJ taking 
the Brillouin number as our central value and the average difference to the Wilson number 
(0.006) as a uniform estimate of the systematic error. Adding all errors in quadrature yields 



M 0++ 

" "ccc 



2.68(2) GeV 
3.71(2) GeV 
4.77(3) GeV 



(18) 
(19) 
(20) 



and we emphasize that these errors do not include the effect of the (missing) limits a — > and 
M^ ea — y 134.8 MeV. Nevertheless, (18) is consistent with experiment, albeit with a large error. 



From these numbers it appears that the actual splitting is very close to equidistant, a notion 
which is also conveyed by Fig.[9| 

We refrain from comparing our numbers to similar results in the recent literature on charm 
physics on the lattice 22 , 24 - 34 . We rather like to add that what is really called for, in our 



opinion, is a complete study with a reasonable assessment of all systematics involved, that is 
with the continuum limit taken, with an interpolation or extrapolation to the physical values 
of M n and Mk in the sea, and with an extrapolation to infinite box volume. 
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6 Summary 



The goal of this work was to test whether a significant difference in meson and baryon dispersion 
relations is seen, depending on whether such a composite state is built from Wilson or Brillouin 
fermions. The main result is that for standard lattice spacings (a -1 = 2 — 3 GeV) this is not the 
case if all quarks are at most as heavy as the physical strange quark, but significant differences 
become visible if one or several quarks are in the range of the physical charm quark mass - see 
Fig.[6]for mesons and Fig.[7]for baryons. 

The Brillouin operator as proposed in 12 can be seen as a low-cost approximation to the 
concept of "perfect fermions" (35 38 . A recent development in the field of staggered fermions 



is to add a mixture of taste-S,V,T,A,P mass terms 39-43 such that the resulting action would 
have only one or two species in the continuum, and an eigenvalue spectrum similar to the 
near-Ginsparg- Wilson spectrum of the Brillouin action (see e.g. Fig. 22 of Ref. [l~2]). 

In summary, we reach the conclusion that the added expense, in terms of CPU time, that 
the Brillouin action entails is hardly justified if one is only interested in light quark spectroscopy 
(this may be different for structure functions). On the other hand, as soon as charm quarks are 
involved, the Brillouin action leads to a massive reduction of cut-off effects already in purely 
spectroscopic quantities. In particular our Fig. [7] shows that the standard lore that aM should 
not exceed 1 need not be true with Brillouin fermions; in this figure we see no cut-off effects 
in the range (aE) 2 = 3.0 — 3.5. All together, it seems our choice to use the Brillouin action to 
determine the quark mass ratio m c /m s in Ref. 44 was justified, and we hope that this augers 
well for the accuracy of our results (19 20). 
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1/2 


-1/2 


c = 2 


s = 


^++ 

' — 'cc 


e(c T C 75 M)c 


1/2 


+1/2 






1 — 'cc 


e(c T C 7 5(i)c 


1/2 


-1/2 


s = 1 


Q + 

cc 


e(c T C 75 s)c 









Table 2: Interpolating fields ofSU(4) 20' -plet ("octet-type") baryons with spin 1/2. Throughout 
the color indices are suppressed, i.e. e{x T C^f 5 y)z is to be read as € a bc( x a^l5yb) z c- 



Appendix: Charmed baryon interpolating fields 



To avoid any confusion to which states our numbers (jl8j|20J) would refer to, we give a list of the 
simplest baryon interpolators with charm and/or strangeness. We follow the naming convention 



of PDG 18 and group the operators according to their transformation properties under SU(4) 
in flavor space. Throughout, C denotes the charge conjugation matrix, the transposition sign 
refers to spinor, and color indices are implicit, as described in the table captions. 

Overall, the states separate into a 20'-plet of spin 1/2 states, a 20-plet of spin 3/2 states, 
and a 4-plet under SU(4). These states are listed in Tab.|2j Tab.[3j and Tab.[IJ respectively. 

The 20'-plet decomposes into the standard c = ground floor which transforms as an 8 
under SU(3) [lines 1-8 in Tab. [2], the c = 1 first floor which decomposes into a 6 [lines 9-14] 
and a 3 [lines 15-17], and the c = 2 second floor which transforms as a 3 under SU(3) [lines 
18-20]. Regarding the c = 1 level, it is worth noticing that the states of the 6 are symmetric 
under interchange of the two non-charmed quarks, whereas the states of the 3 are antisymmetric 
under this interchange. Here we adopt the rule that in the 20'-plet the diquark {x T C r )^y) is 
antisymmetric under the interchange x +>■ y. 

The structure of the 20-plet is somewhat simpler, since each fixed-c floor has a unique trans- 
formation pattern under SU(3). It contains the standard c = ground floor which transforms as 
a 10 under SU(3) [lines 1-10 in Tab.pl, the c = 1 first floor which transforms as a 6 [lines 11-16], 
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charm strange 


baryon interpolating field / I z 


c = s = 


A++ e(ti r C7 M «)« 3/2 +3/2 
A+ ^e{2(w T C 7 ^)w + (m t C 7/ ,m)4 3/2 +1/2 
A ^e{2(d T C 7At u)d+ (d T C 7M rf)M} 3/2 -1/2 
A- t{d T C lfl d)d 3/2 -3/2 


s = 1 


£*+ ^e{2(w r C 7M s)w + w r C 7M w)s} 1 +1 
£*° ^ e {( M T C7^)s + (d T C 7/ ,s)u + (s T C7 M u)4 1 
£*~ X ^e{2(d T Cj IM s)d + d T Cj t ,d)s} 1 -1 


s = 2 


S*° j=e{2(s T C lfl u)s + (s T CVM 1/2 +1/2 
S*~ l f {2(/C 7 ,<l) S + (/Cv)(i} 1/2 -1/2 


s = 3 


0~ e{s T C-f fl s)s 


c = 1 s = 


£*++ ^e{2(« T C 7M c) M + M T C 7AtM )c} 1 +1 
K + ^e{(w T C 7/ ^)c + (d T C lfl c)u + (c T C 7 ^)rf} 1 
£*° V> ^e{2(rf T C 7/ ,c)rf + cFCj^c} 1 -1 


s = 1 


n*+ ^e{(« T CV0c+ (s t C 7m c)m + (c t C 7/ ,m)s} 1/2 +1/2 
^e{(d T C 7MS )c + (s T C 7 ^c)d + (c T C lfl d)s} 1/2 -1/2 


s = 2 


Of ^e{2( S r C 7/ ,c) S + (s r CV)c} 


c = 2 s = 
s = 1 


H* c ++ ^e{2(c 1 C 1 ,u)c+(c 1 C ltl c)u} 1/2 +1/2 
e£" ^e{2( C T CV)c+(c T C 7/ ,c)4 1/2 -1/2 
02+ ^e{2(c r C 7 ^)c+(c r CVM 


c = 3 s = 


0++ 6(c t C 7m c) C 



Table 3: Interpolating fields of SU(4) 20-plet ( u decuplet-type" ) baryons with spin 3/2. 
Throughout the color indices are suppressed, i.e. e(x T Cj^y)z is to be read as ea^x^Cj^y^Zc- 



charm strange 


baryon 


interpolating field 


/ 




c = s = 1 


A'° 


^e{{u r C lb d)s + (s r C lb u)d + (d r C lb s)u} 








c= 1 s = 


K + 


^e{(u r C lb d)c+ (c r C^u)d+ (d r C lb c)u} 








s = 1 


' — 'c 
' — 'c 


^e{(u T C l5 s)c + (c r C l5 u)s + (s t C 75 c)m} 
^e{(tFC7 75 s)c+ (c T C l5 d)s + (s T C 75 c)4 


1/2 
1/2 


+ 1/2 
-1/2 



Table 4: Interpolating fields of SU(4) 4-plet ( u new-type" ) baryons. Throughout the color 
indices are suppressed, i.e. e{x T C^^y)z is to be read as e a 6 c (a;^ C^yi,) z c . 



the c = 2 second floor which transforms as a 3 [lines 17-19], and the c = 3 one-point summit 
of the pyramid [line 20]. Here we adopt the rule that in the 20-plet the diquark (x T C r y^ l y) is 
symmetric under the interchange x -<->■ y. 

The 4-plet decomposes into a c = ground floor which is an SU(3) singlet [line 1 in Tab. [I], 
and a c = 1 first floor which transforms as a 3 [lines 2-4]. In the former case the construction 
is based on the requirement that S° oc (us)d + (ds)u, A = Ag oc 2(ud)s + (us)d — (ds)u, 
and A'° = Aq oc (ud)s + (su)d + (ds)u would be mutually orthogonal. In the latter case the 
interpolator is antisymmetric under the interchange of the two non-charmed quarks, if we adopt 
the rule that the diquark {x T C r )^,y) is antisymmetric under the interchange x ■<-»• y. 
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